Frustration-Driven Successive Metamagnetic Transitions in TbB 4 
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Resonant magnetic x-ray diffraction experiments on the Shastry-Sutherland lattice TbB4 were 
carried out under strong pulsed magnetic fields up to 30 T. TbEU exhibits a multi-step magnetization 
process above 16 T when magnetic fields are applied along the c-axis. We examined the intensity of 
the 010 magnetic reflection as a function of magnetic field and found that the magnetization plateau 
phases are accompanied by large XY components of magnetic moments, in contrast to normal 
fractional magnetization plateau phases. The magnetization was calculated using a simple spin 
model deduced from the above result. Finally we propose that frustration is the key to understanding 
the observed multi-step magnetization. 

PACS numbers: 75.25.+z,75.30.Kz,78.70.Ck 



When one applies pressure to a substance, the volume 
will be reduced at all times depending on its compress- 
ibility. In contrast to such natural behavior, several mag- 
netic materials exhibit magnetization plateaus, where the 
magnetization keeps a constant value at a finite range 
of magnetic fields even below the saturation field. This 
"zero compressibility" is simply ascribed to the fact that 
S z is a good quantum number, while a variety of mech- 
anisms give rise to plateaus in the magnetization curve. 

Let us consider only plateau phases in ordered states 
here. The easy-axis magnetization of an Ising antiferro- 
magnet is most comprehensible. An example is the one- 
third magnetization plateau of CoCl2-2H20 [l,]. Compet- 
ing exchange interactions stabilize the fff structure just 
below the saturation field. The multi-step magnetization 
curves of PrCo2Si2 and CeSb Q are also classified 
into this category. The plateau phases are described by 
a complex sequence of up- and down-spins. 

Magnetization plateaus are also observed in Heisenberg 
spins. In geometrically frustrated Heisenberg antiferro- 
magnets, one-third (ft J,) and one-half (ffff) plateaus 
are sometimes observed in triangular- and tetrahedon- 
based lattices, respectively. Although the field range of 
the plateau phase is zero for these lattices if only the 
nearest- neighbor interaction is taken into account, ad- 
ditional mechanisms, such as a four-spin exchange cou- 
pling, quantum fluctuations, single-ion anisotropy and 
magnetoelastic couplings, stabilize the plateau phase in 
real magnets 0,0, S0- A crucial point is that all plateau 
phases have a collinear magnetic structure. No ordered 
magnetic moment perpendicular to the magnetic field ex- 
ists in the plateau phases. 

Recently, a multi-step magnetization curve was ob- 
served in TbB4. Intriguingly, several experimental evi- 
dences indicate that TbB4 resembles none of the above- 
mentioned magnets. In this letter, we report on reso- 
nant x-ray diffraction (RXD) experiments on TbB 4 up 



to 30 T. The field dependence of the magnetic intensity 
provided decisive information on the magnetic structure 
of the field-induced phases. From this result, we show 
that TbB 4 indeed belongs to a rare class of magnets, 
characterized by ordered perpendicular components in 
the plateau phases. We also show that frustration plays 
a vital role in the multi-step magnetization process on 
the basis of the calculation of magnetization curves. 

The rare-earth tetra-boride i?B 4 crystallizes into a 
tetragonal structure (space group: PA/mbm). The net- 
work of the rare-earth ions is equivalent to the frustrated 
Shastry-Sutherland (SS) lattice [§]. Several pieces of 
work have thus been carried out from a viewpoint of frus- 
tration between rare-earth quadrupole moments 0, ITo| . 
Somewhat confusingly, in the strict sense, classical con- 
tinuous spins on the SS lattice are not frustrated, because 
the ground state is a unique helical or Neel structure and 
has no degeneracy [8(. However, for helical structures, 
none of the antiferromagnetic bonds are antiparallel and 
thus local pair exchange-energy is not minimized. The 
spins are not satisfied and still feel "frustration" . We 
will revisit this point later. 

In TbB 4 , there are two transition points Tni=44 K and 
Tn2=22 K. The total angular momentum of the Tb ion is 
large (J=6) and hence the Tb spin is approximated by a 
classical spin. The magnetic structure between Tni and 
Tn2 revealed by neutron powder diffraction experiments 
is shown in Fig. [TJa) [llj. The magnetic moments are 
XY-type and are confined in the basal afe-plane. Addi- 
tional in-plane anisotropy orients the moments parallel 
to the diagonal line. The nearest-neighbor interaction Ji 
is antiferromagnetic. The next-nearest-neighbor interac- 
tion Ji is probably antiferromagnetic, and J\ and J2 form 
a frustrated SS lattice. In i?B 4 , the further-neighbor anti- 
ferromagnetic interaction J3 is taken into account, which 
stabilizes the observed magnetic structure. In this struc- 
ture, the Ji bond is not satisfied. Below Tn2, the mag- 
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FIG. 1: (a) Magnetic structure of Tbl>4 between Tni and Tn2- 
A large square is the unit cell. Major interactions ( J\ : broken 
line, J2 : solid line, J3 : dot-dashed line) are shown. Sample 
coordination (a, b, c, x, y and z) and incident and scattered x- 
rays (fci and fcf) are also shown, (b) Magnetization process of 
TbB 4 taken from Ref. QJ]. 



netic moments rotate in the afr-plane and a collinear-like 
antifcrromagnetic structure is realized. This transition 
is considered to be a ferro-quadrupole order [13] and is 
destroyed by application of magnetic fields. 

The magnetization process of TbB 4 parallel to the c- 
axis shown in Fig.[ljb) exhibits multi-step metamagnetic 
behavior between 16 and 28 T The number of the 
steps is reported to be nine, and the one-half magne- 
tization plateau is most prominent. When a magnetic 
field is applied parallel to the [100] and [110] directions, 
a single metamagentic transition is observed. Magneti- 
zation measurements of diluted samples indicated that 
the multi-step magnetization is inherent to the phase be- 
tween Tni and Tn2 



14| . The metamagnetic behavior 



suggests that the essential nature of the transitions is 
level-crossing between low- lying multiplets [ToT ]. How- 
ever, the origin of the multiple successive transitions is 
not easily understood. Magnetic structures in the high- 
field plateau phases provide useful information, thus we 
conducted an x-ray diffraction experiment under strong 
pulsed magnetic fields, which has been recently devel- 
oped utilizin g b rilliant third-generation synchrotron x- 

rays HE [13, Hi. Althou g h 

magnetic x-ray scattering 
cross-sections are very small, we utilized resonant en- 
hancement of magnetic scattering [lij] . 

Synchrotron x-ray diffraction experiments were carried 
out at undulator beamline BL22XU at SPring-8. The 
x-ray energy was tuned to the Tb Tni absorption edge 
(7.514 keV). Magnetic fields up to 30 T were generated 
using a small pulsed magnet (20 mm in outer diameter 
and 24 mm in length). The pulse duration was about 
0.6 ms. The magnet was a split-pair, and two windows 
were prepared for incoming and outgoing x-rays. The 
magnet and a sample were attached to the 100 K and 
10 K stages of a conventional closed cycle refrigerator, 



respectively. Details of experimental setup are described 
in Ref. [l6|. Single crystals of TbB 4 were grown by the 
floating zone method. The sample was a thick plate of 
dimensions 1 mmxl mmx0.5 mm with a polished (010) 
surface. The FWHM of the u profile of the 020 reflection 
was about 0.06°. 

The scattering plane was horizontal. The b* axis lay 
in the scattering plane and the c-axis was perpendicular 
to the scattering plane. We observed the forbidden 010 
reflection, which is a magnetic Bragg point as demon- 
strated in the neutron experiment [ill ] . The azimuthal 
angle \&, which is the rotation about the scattering vec- 
tor, is a key parameter in RXD. We define <F=0 when 
the c-axis is perpendicular to the scattering plane. 'F 
was kept close to zero throughout this experiment. Mag- 
netic fields were applied parallel to the c-axis. We used 
a high-counting-rate detector and measured diffraction 
intensity as a function of time with a typical resolution 
of 10 /is. The field strength was also monitored as a 
function of time, and we obtained the field dependence 
of the diffraction intensity from these two series of data. 
Polarization of scattered x-rays were separated using a 
PG (006) crystal analyzer. In order to prevent the sam- 
ple from being heated by x-ray irradiation, we inserted 
a chopper into incident x-rays. The window width was 
about 2 ms and the duty cycle was about 0.01. 

In Fig. [2ja), we show the peak intensity of the 010 re- 
flection as a function of photon energy below Tn2 and 
above Tni- The fluorescence spectrum is also shown. 
The intensity of the 010 reflection exhibits huge enhance- 
ment at the L3 main edge, implying that the resonance 
is ascribed to electric dipole (El) transition. Subse- 
quent measurements were carried out at the peak energy 
(7.512 keV). The temperature dependence of the inte- 
grated intensity of the 010 reflection is shown in Fig.^b). 
The intensity appears below Tni and is constant below 
Tn2- The effect of beam heating is estimated to be at 
most 2 K from the differences between the reported and 
observed transition temperatures. 

Here we briefly mention the resonant x-ray scattering 
amplitude [2l| of the Tb ion, which is proportional to 



/ rcs oc iCi (e| x 6i) • TTl + C 2 e\ O 6i 



(1) 



where m is the magnetic moment and Cj and ef are 
the polarization vectors of the incident and scattered x- 
rays, respectively. The symmetric second-rank tensor 
O describes anisotropy of the Tb 5d orbital caused by 
anisotropic crystal environment or quadrupole order and 
is represented by a linear combination of the five elements 
O xy , O yz , O zx , O22 and O2o- The first term in eq.{T|) is 
magnetic scattering and the last term causes anisotropy 
of the tensor of susceptibility (ATS) scattering. Thor- 
ough measurements on temperature and azimuthal angle 
dependence with polarization analysis were performed by 
TM [2(j. The results were consistent with the magnetic 
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FIG. 2: (color online) (a) Peak intensity of the 010 reflection 
as a function of photon energy below Tn2 and above Tni. 
Huge enhancement is evident at the L3 edge of the fluores- 
cence spectrum. Small peaks observed at 50 K are multiple 
scattering, (b) Intensity of the 010 reflection as a function of 
temperature. Arrows indicate Tni=44 K and Tn2=22 K. The 
sample temperature is higher than the thermometer by 2 K 
because of beam heating, (c) Field dependence of the peak 
intensity of the 010 reflection at 8 K. The arrows indicate the 
onset of the high-field phases (16.4 T), the lower and upper 
critical fields of the one-half plateau phase (20.5 T and 24.7 T) 
and the saturation field (28 T) at 4 K. (d) Field dependence 
of the 010 reflection at 8~10 K. Polarization of scattered x- 
rays are separated. Red and blue (magenta and cyan) circles 
are up- and down-sweeps of the n-a (n-n) channel. 



structure obtained in the neutron experiment and 
it was confirmed that OfcO reflections (k : odd integer) at 
^=0 are magnetic scattering. 

The field dependence of the peak intensity of the 010 
reflection at 8 K is shown in Fig. [^c) . The data were av- 
eraged over 12 field-scans. The intensity is constant up 
to 16 T and gradually decreases above 16 T as the field 
is increased. An important outcome is that a consider- 
able amount of the intensity remains in the magnetization 
plateau phases above 16 T. From this result, we under- 
stand that the plateau phases include (0,1,0) modulated 
magnetic moments and/or anisotropic 5d orbitals. 

The second important outcome is obtained from polar- 
ization dependence of scattered x-rays. States of linear 
polarization are labeled by a and 7r when they are normal 
and parallel to the scattering plane, respectively. The 
incident polarization was 7r. The scattering amplitudes 
f rcs {9) of tttt' and ixa' channels for each magnetic and 
ATS component at \P=0 are listed in Table [H where 9 is 
the Bragg angle. The structure factor f rcs (9) e l2vKRi 
at (0,fc,0) is proportional to either if res (9) sin 2irk5 or 
f ICS (9) cos 27rfc<5, where the fractional coordinate of the 
Tb site S ~ 0.317. The structure factor below 16 T is 
proportional to m x cos#sin27rfc(5. In Fig. HJd), field de- 
pendence of 7T7r' and na' channels of the 010 reflection 



TABLE I: Scattering amplitude of each component of mag 
netic and ATS scattering at 9=0 for titt' and txcj' channels. 



are shown. The data were averaged over 30 field-scans. 
It is obvious that the scattered x-rays are dominated by 
7rer' channel at all field ranges. Hence it is inferred from 
Table[T]that the scattering intensity of the 010 reflection 
arises from m x +m y or O yz +O zx . We also observed the 
050 reflection (not shown) without the analyzer. In spite 
of large difference in 9 (6.6° and 35.2°), the field depen- 
dence of the 050 reflection is very similar to that of the 
010 reflection. This fact indicates that the intensities of 
both the zero- and high-field phases are proportional to 
cos 9 sin 27rfc<5 and that m z (oc sin 29) does not exist in the 
high-field phases. In addition, recent neutron diffraction 
experiments of TbB4 under pulsed magnetic fields have 
shown that the field dependence of the 010 reflection in 
neutron diffraction is also similar to that of x-ray diffrac- 
tion [221 ] . Neutrons observe m x , m y and m z . Hence pure 
Oyz+Ozx order is inconsistent with the neutron result. 
m x or m y must exists. Actually all data are well inter- 
preted by assuming that the major order parameter in 
the high-field phases is m x , magnetic moments perpen- 
dicular to magnetic fields. 

This is a rather surprising result, because all magnetic 
moments in plateau phases are parallel or antiparallcl to 
the magnetic field as mentioned in the introduction. A 
canted magnetic structure, usually derived from the co- 
existence of parallel and perpendicular components, re- 
sults in a continuous magnetization curve unlike the ob- 
served step-like magnetization process. Accordingly, we 
submit an idea that in TbB4 the magnetic structures in 
the plateau phases consist of very hard XY-type magnetic 
moments and Ising-like magnetic moments. 

In order to interpret this novel magnetic structure, we 
propose the following spin model. The low-lying mul- 
tiplets are approximated by the ground-state XY-spin S 
with the in-plane anisotropy E and the excited-state Ising 
spin s with the energy gap G. The spin Hamiltonian is 



h = E -M E ( W + s * s j) + E ^ 

ra=l,2,3 V i,j k,l 



where the spins are unit vectors, i and j sites are occupied 
by XY spins and k and I sites are occupied by Ising spins. 

Magnetization curves were calculated for finite-size 
cells with periodic boundary conditions. We consider 4 
spins in the crystallographic unit cell (lxl cell). For 
simplicity, we initially consider a four-clock model by 
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FIG. 3: (a) Calculated magnetization curve for E — > oo. The 
magnetic structure of the one-half plateau phase is shown in 
the inset. denotes a spin parallel to the 2-axis. (b) Possible 
magnetic structure of the one-third plateau phases. 



taking the limit E — > oo. For example (see Fig. Efa)), 
allowed states for spins 1 and 3 are / ( — 2 ' — 2 >0)> 
/ (^,^,0), 0(0,0,1) and ®(0,0,-l). In this case, we 
calculated energies of all 4 4 =256 states and found ground 
states. For rather wide parameter sets, this model has 
the same ground state as observed. A calculated magne- 
tization curve for J\ — 2, J2 — 1, J3 = 2 and G = 1 is 
shown in Fig. (Ha). The one- half magnetization plateau 
is clearly reproduced and the corresponding magnetic 
structure is shown in the inset. Similar results were ob- 
tained with other sets of parameters and different cell 
sizes (2x2 and 1x3 cells). 

It is interesting to note that a number of states that 
have various magnetizations are degenerate at the criti- 
cal fields H c i/2 and H c \. This is in fact an ordinary phe- 
nomenon observed in metamagnetic transitions. How- 
ever, if this degeneracy is lifted, the multi-step magneti- 
zation curve of TbB4 can be reproduced. Three mech- 
anisms can lift the degeneracy: (i) difference in magni- 
tude between S and s, (ii) interactions beyond J3, and 
(iii) frustration. While the first two mechanisms cause 
symmetric splits at H cl / 2 and H c \, the observed mag- 
netization curve is asymmetric. For instance, the one- 
third plateau is clearly visible, whereas the two-third 
plateau is indistinguishable, as seen in Figfljb). It is 
evident that the last mechanism, frustration, plays an 
essential role. A possible magnetic structure of the one- 
third plateau phase is shown in Fig. [3^b). Here we set 
the in-plane anisotropy E finite so that the moments can 
rotate in the a6-plane. The Ising sites behave like de- 
fects. In frustrated magnets, as pointed out in the in- 
troduction, the local pair exchange-energy is not mini- 
mized. Hence, when a defect is introduced, the magnetic 
structure around the defect relaxes and gains exchange 
energy. This does not happen in magnetic lattices with- 
out frustration, because all exchange energies have been 
minimized. 

For magnetic systems with short-range interactions, it 
is expected that small commensurate unit cells are more 
stable than large or incommensurate cells. Thus it is 



plausible that a finite number of selected magnetization 
steps appear in the magnetization curve through the frus- 
tration mechanism. This mechanism is active only when 
the density of the Ising sites is less than half (up to J3). 
Therefore the degeneracy is lifted only at H cl / 2 - The 
diffraction intensity of magnetic structures based on the 
model shown in Fig. [3fa) is one-quarter of the intensity 
at T. The observed intensity is, however, about 40% of 
that at T, implying presence of other order parameters, 
such as O20 and O yz . 

In conclusion, we have investigated the multi-step mag- 
netization process of TbB4 using RXD. The magnetic 
intensity observed in the plateau phases illustrates that 
large XY-components of the magnetic moments coexist 
with aligned Z-components. A simple spin model repro- 
duces the one-half magnetization plateau, and it is found 
that the frustrated nature of the SS lattice lifts the de- 
generacy at the critical fields and creates the observed 
multi-step metamagnetic behavior. Paradoxically, frus- 
tration is not only an origin of degenerate states, but also 
breaks degeneracy. 

This work was supported by Grant-in- Aid for Scientific 
Research on Priority Areas "High Field Spin Science in 
100T" (No.451) from the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT). 



[1] 

[2. 
[3; 
[4 

[5 
[6 

[7] 
[8] 

[9 
[10 

[11 

[12 

[13 
[14 
[15 
[16 

[17 
[18 

[19 
[20 
[21 

[22] 



Phys. Soc. Jpn. 74, 2169 (2005). 
Phys. Soc. Jpn. 74, 2434 (2005). 
, J. Phys. Soc. Jpn. 76, 015001 



H. Kobayashi and T. Haseda, J. Phys. Soc. Jpn. 19, 765 
(1964). 

T. Shigeoka et al., J. Phys. Soc. Jpn. 58, 394 (1989). 
J. Rossat-Mignod et al., Phys. Rev. B 16, 440 (1977). 
T. Sakakibara and M. Date, J. Phys. Soc. Jpn. 53, 3599 
(1984). 

L. E. Svistov et al., Phys. Rev. B. 67, 094434 (2003). 
H. Kitazawa et al, Physica B 259-261, 890 (1999). 
H. Ueda et al., Phys. Rev. Lett. 94, 047202 (2005). 
B. S. Shastry and B. Sutherland, Physica B 108, 1069 
(1981). 

R. Watanuki et al., J. 
D. Okuyama et al., J. 
T. Matsumura et al. 
(2007). 

T. Fujita et al., Abstr. Meet. Phys. Soc. Jpn. (63rd Annu. 

Meet., 2008), Part 3, 23aPS-86 [in Japanese] . 

S. Yoshii et al., J. Magn. Magn. Mater. 310, 1282 (2007). 

F. Iga et al., private communications . 

Z. A. Kazei et al., JETP Letters 84, 441 (2006). 

Y. H. Matsuda et al., J. Phys. Soc. Jpn. 75, 024710 

(2006). 

Y. Narumi et al., J. Synch. Rad. 13, 271 (2006). 
P. Frings et al., Rev. Sci. Instr. 77, 063903 (2006). 
D. Gibbs et al., Phys. Rev. Lett. 61, 1241 (1988). 
T. Matsumura et al., private communications . 
M. Blume, Resonant Anomalous X-Ray Scattering (El- 
sevier Science B.V., Amsterdam, 1994). 
S. Yoshii et al., private communications . 



